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Summary 

A theory to explain the initial stages of unsteady separation has been proposed by Van 
Dommelen Sc Cowley (1989). In the present paper, this theory is verified for the separation 
process that occurs at the equatorial plane of a sphere or a spheroid which is impulsively 
spun around an axis of symmetry. A Lagrangian numerical scheme is developed which 
gives results in good agreement with Eulerian computations, but which is significantly 
more accurate. This increased accuracy, and a simpler structure to the solution, also 
allows verification of the Eulerian structure, including the presence of logarithmic terms. 
Further, while the Eulerian computations broke down at the first occurrence of separation, 
it is found that the Lagrangian computation can be continued. It is argued that this 
separated solution does provide useful insight into the further evolution of the separated 
flow. A remarkable conclusion is that an unseparated vorticity layer at the wall, a familiar 
feature in unsteady separation processes, disappears in finite time. 

1. Introduction 

In part 1, Van Dommelen & Cowley (1989) proposed a Lagrangian description for 
unsteady separation under a wide range of conditions. This theory has been verified in a 
number of two-dimensional unsteady boundary-layer computations, the first being the one 
by Van Dommelen Sc Shen (1980,1982) for a circular cylinder which is impulsively set into 
motion in the direction normal to its axis. Yet it seems somewhat unsatisfactory that the 
verification of basic aspects of the theory should depend solely on two-dimensional unsteady 
computations, since by necessity their resolution is much lower and their convergence 
questions more complex than one-dimensional computations. This motivated the present 
examination of the separation process which occurs at the equatorial plane of a sphere 
which is impulsively rotated. 
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This separation process was first described by Banks &: Zaturska (1979), who proposed 
an analytical description in the form of a power series in time. However, Simpson Sc 
Stewartson (1982) argued that there would also be higher order logarithmic terms in the 
expansion. (Interestingly, Banks & Zaturska (1981) discovered the presence of logarithmic 
terms in a flow previously studied by Bodonyi &; Stewartson (1977), cf. part 1). The 
physical separation process was explained in Lagrangian terms by Van Dommelen (1981), 
who proposed that separation would cause the boundary layer to divide into an unseparated 
vortex layer at the wall and a separated one above it. His proposals agree with the Eulerian 
expansions of Simpson & Stewartson (1982). While all these studies were only concerned 
with the flow in the equatorial plane, numerical solutions for the complete boundary-layer 
flow were given by Dennis & Ingham (1979) and Van Dommelen (1986). Navier-Stokes 
solutions were presented by Dennis & Duck (1988). 

Except as an example of a basic unsteady separation process, the spinning sphere is 
also of interest because of its relation to problems in geology and meteorology. It further 
turns out that the results apply equally well to axially symmetric bodies of general shape, 
provided that the body is symmetric about the equatorial plane (cf. part 1 and Banks & 
Zaturska 1979). 

The structure proposed by Van Dommelen & Cowley (1989) is verified in sections 
3 and 4. For this purpose, in section 2 a Lagrangian numerical procedure is described. 
This procedure is found to be more accurate than the Eulerian Crank-Nicolson scheme of 
Banks & Zaturska (1979) and the Eulerian box scheme of Simpson & Stewartson (1982). 
It allows a precise verification of the presence of logarithmic terms in the expansions 
as proposed by Simpson & Stewartson (1982) (section 5). Unlike the Eulerian schemes, 
the Lagrangian computation can be continued beyond the first occurrence of separation 
without apparent difficulties. The physical meaning of such a solution is not immediately 
clear, since interaction effects invalidate the equations in the immediate vicinity of the 
equatorial plane. Yet, in section 6 we will argue that our solution can be extrapolated 
away from the equatorial plane to where the governing equations are still correct. Thus 
our solution may provide an interesting first glimpse of the continued evolution of an 
unsteady separated flow. One remarkable result is that the unseparated vortex layer at 
the wall disappears quickly. In section 7 an analytical description for this process is derived 
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which is in good agreement with the numerical data. 

2. Lagrangian Formulation And Numerical Method 

The flow about the spinning sphere is most easily described in a polar coordinate 
system where x is the polar angle measured from the axis of rotation and y the distance 
from the wall, scaled to eliminate the coefficient of viscosity. The corresponding velocity 
components are u and v, and w denotes the velocity in the azimuthal direction. The radius 
of the sphere and the inverse of the angular velocity are scaled to unit values. 

Van Dommelen & Cowley (1989) further introduce Lagrangian coordinates £ and 77 
attached to the fluid. By definition, these are taken equal to the values of x and y at the 
time, t = 0 , that the sphere initially starts to spin. In terms of £ and 77 , the boundary-layer 
momentum equations, obtained by a series expansion in f at the equator, are: 

x ,{ = u ,« . (la) 

w = x^w i ,„ +x t( x )(v w t1} , (16) 

^ - X ,« X *™K { - x if w 2 , (lc) 

where the subscript comma denotes differentiation with respect to the subsequent sub- 
scripts, the dot denotes differentiation with respect to time, and x^, u^, and w are eval- 
uated at the meridional plane £ — 7 r/ 2 , so that they depend on 77 and t only. The first 
equation, (la), is the ^-derivative of u — x which defines the polar flow velocity as the 
time derivative of the polar position, (16) is the azimuthal momentum equation as given 
by Van Dommelen & Cowley (1989), while (lc) is the ^-derivative of the polar momentum 
equation. 

The initial and boundary conditions are: 

= 1 » u ,(( 7 7 ) 0 ) = 0, (1 d,e) 

w(t/,0) =0 if V ± 0 , ( 1 /) 

x ,f(M) = 1 , u ,{(0> t) = 0 , w(0,t) = 1 , (I g,h,i) 

x e (oo, t ) = l, u ( (oo,t)= 0 , w(oo,t) = 0 . (1 
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The advantage of Lagrangian coordinates arises from the fact that the particle distance 
from the wall, y, occurs only in the continuity equation (cf. part 1): 


r h 

0 t) ’ 


( 2 ) 


This equation needs to be integrated only at times at which results are desired; it does 
not affect the numerical solution of (la) through (1/). It also turns out that the continuity 
equation is the first one to become singular, so that numerical difficulties do not arise in 
the integration of (la) through (1/). 

The present numerical integration follows the general lines of the procedure of Van 
Dommelen & Shen (1980). For example, to achieve an effective distribution of mesh points 
across the boundary layer, an arctangent mapping was used, and the singularity at the 
impulsive start was eliminated by a further coordinate transformation. The Jacobian (2) 
was integrated using quadratic interpolation for . 

The Lagrangian equations (1) were discretisized by means of Crank-Nicolson central 
finite differences. The resulting implicit finite difference equations were solved iteratively 
for xt, w, and u^ respectively by means of the tridiagonal algorithm. The iterations 
were continued until the error in the finite difference equivalent of (la) through (lc) was 
less than 3 10~ 8 ; this avoids the possible problems of a termination criterion based on 
the difference between iterates. To eliminate possible round-off errors, 17 digit numerical 
precision was used throughout. 

Thus the major source of numerical inaccuracy should be the truncation error. To 
take account of this error, computations were performed at the four meshes listed in table 
1, which compare favorably to the meshes used in the Eulerian computations. 

In contrast to the Eulerian calculations, no Richardson extrapolation was used. But 
if desired, it may be noted that the Lagrangian solution should provide a much better 
basis for repeated Richardson extrapolation than the Eulerian schemes. The reason is the 
smoothness of the Lagrangian solution discussed in the next section. 


3. Separation Structure 

In solving the problem specified by equations (1), (2) in Enfrrian coordinates. Banks 
Zaturska (1979) discovered that the boundary-layer thickness becomes infinite at some 
finite time ts- According to Sears & Telionis (1975), such singularities in a classical 
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boundary-layer solution indicate separation. The singularity should be understood to 
mean (hat the loral motion away from the wall becomes too strong to be described with 
boundary-layer scalings. As an example, Elliott, Cowley &: Smith (1983) show that the 
interactive stage of unsteady two-dimensional separation occurs at a boundary-layer thick- 
ness 0(Re ~ 5 / 11 ), rather than the classical 0(Re~ 1 ^ 2 ). 

The separation processes proposed by Van Dommelen & Cowley (1989) axe character- 
ized by non-singular solutions x,(, u^, and w to the Lagrangian boundary-layer equations 
(1). Singular behaviour should occur only in the continuity equation (2), caused by vanish- 
ing of x t ( inside the boundary layer. Such behavior results in infinite values of the particle 
position y, leading to the infinite boundary-layer thickness observed in the Eulerian com- 
putations. 

The numerical results do show that the Lagrangian x^, u^ and w profiles remain 
regular, and that x,( becomes zero. The first zero occurs at a point S located at tjs — 
0.97188 at time ts = 4.575632. This time is in excellent agreement with the time that 
the Eulerian solution becomes singular, 4.5758 according to Banks & Zaturska (1979) or 
4.57446 according to Simpson & Stewartson (1982). (Values of the computational quantity 
used by Simpson Sz Stewartson to find the separation time are listed in table 2). Yet unlike 
the Eulerian separation profiles, the Lagrangian profiles of figure 1 do not show any sign 
of singular behaviour. 

The vanishing of x^ leads to a singular solution to the continuity equation (2) for the 
y— position of the particles. To find the structure of this singularity, x t ^ can be expanded 
in a finite Taylor series expansion around the point S: 

x ,t = 5 x ,tvv(v'A s )6t) 2 +u i (-n,t')Si . (3a) 

The precise requirement for this series to be valid is that the appearing derivatives are well- 
defined, i.e. continuous near the point S. The computed x^,,,, and u t { profiles at time ts, 
figure 1, show no sign of singular behaviour at point S, and prove highly accurate according 
to comparisons for varying mesh size. As an example, table 3 lists the convergence of the 
values of x.,( vv and u^ at the point S, along with their ^-derivatives. 

The Taylor series expansion (3a) may be substituted into the continuity equation (2) 
to find the vertical position y of the particles: 

Y = M *y = | (axetan + £) + 4 1 ) , ( 3b ) 
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V = 0(1) • 


f3 = (-2x^u { )= , 


(3c, d) 


Here, the omission of the subscript comma indicates that the value of the derivative at the 
point S is meant. This paper follows the Eulerian definitions of the constants closely; in 
terms of the notation used in part 1, /? corresponds to and /3/x^ vv to /? 0 /? 2 . Our 

value for the constant /? = 0.71387, table 3, is in good agreement with the 0.71 of Banks 
Zaturska (1979) and the 0.712 of Simpson & Stewartson (1982). 

Figure 2 shows Eulerian velocity profiles in terms of the above scaled coordinate Y 
and the scaled Eulerian velocity gradient 


G = -|£t|u )X , 



(4a, 6) 


In these velocity profiles, the particles rj < Tjs are found in a wall layer near Y = 0 (see 
(3b)); similarly the particles tj > rjs axe located in a separating layer near (3Y = 27 r. The 
part 0 < /? Y < 27 t comprises most of the Eulerian velocity profile, yet it corresponds to 
only a small vicinity Stj = 0(|£t|s) of the point S in the Lagrangian profile. 

The last observation implies that for 0 < Y < 2n, a Taylor series expansion for 
the Lagrangian solution is applicable. The Lagrangian coordinate rj may subsequently be 
eliminated in favor of Y by means of (36) to find asymptotic expressions for the Eulerian 
velocity profiles: 


w ~ ws , 

u,x~-|ar 1 A(l-cos(/?Y)) . 


(5a) 

(5b) 


The scaling (36) of the variable Y compensates for the rapid expansion of the region 
of particles near tjs ■ This scaling leads to the apparent thinning of the wall layer and the 
separating layer in the velocity profiles of figure 2: in terms of the original coordinate y, 
these two layers remain of finite thickness. 

To find the separation structure to higher order of approximation, it is more convenient 
to replace the finite Taylor series expansion (3a) in favor of a formal matched asymptotic 
expansion. The proper inner ^-coordinate E near the particle rjs can easily be found by 
applying Van Dyke’s (1975) guiding principles: clearly the reason for the non-uniformity 
in the continuity equation (2) is the vanishing of x^ at the particle 77 s - Removal of this 
non-uniformity requires that the time-dependent term in the Taylor series expansion (3a) is 


6 



retained in the inner region. On the other hand, the matching with the wall and separating 
layers can only be done when the second order term is retained also. Therefore an inner 
coordinate E* will be defined as 


St] = \6i\t— E* , 

x tvv 


(6a) 


which is also consistent with (36). An inner dependent variable AY* can be defined as 


y = 




= ys(t) + 


1 


|$t|M 


?AY* 


( 66 ) 


where ys(t) denotes the vertical distance of the particle t]s from the wall, a distance 
which has been subtracted in order to avoid the appearance of logarithmic terms in the 
Lagrangian inner expansions. 

With the inner scalings known, it is a simple matter to integrate (2) to find the particle 
position: 


AY* = arctan(E*) + |6t|j 


(-0B ln(l +E* 2 ) + 


((3B + C)E* 2 \ 
1 -f E* 2 ) 


+ • • • 5 


(7a) 


^ = , JiSf, 

ox (vv p 


(76, c, d) 


The Eulerian velocity profiles axe found by a Taylor series expansion of the Lagrangian 
profiles, followed by the elimination of the Lagrangian coordinate using (7a): 

G ~ |(1 +cos2AY*) + |6t|M^sin(2AY*)ln(i(l +cos2AY*)) +... , ( 8 a) 


w ~ W5 + |6 t| s fth. tan AY* . (86) 

x (vv 

The above Eulerian results agree with the inner structure of Banks Zaturska (1979) 
as modified by Simpson fc Stewartson (1982). 

4. The Wall And Separating Layers 

The inner solution derived in the previous section can be matched below, at E* — -oo, 
to a wall layer of particles tj < rfc, and above, at E* = oo, to a separating layer of particles 
77 > 775 . In the wall layer, all variables including the particle position y are non-singular 
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since the integral (2) does not involve the singular point tjs. The asymptotic description 
of the wall layer is therefore a Taylor series expansion in time: 


1 

(*,O u ,O w >y)= J2 rf \ 6t \ n ( x t^v),»t\v), win Hvhy (n Hv)) f °r*7<*7s . (9a, b) 


It is not possible to add singular terms of the general form |6t| Q In |6t| to this expansion, 
for substitution into the Lagrangian equations (1) and (2) would lead to inconsistencies. 

The Eulerian asymptotic expansion of the wall layer is found by formal elimination of 
tj in favor of y: 


(u x ,v,w)= ^2 ^!^ t r( u ,*x (n) (y)> v * (n) (y)- w * (n) (y)) for y = o{\) . (iOa,&) 

n=0 


The vertical velocity component v is simply the integral of u )X with respect to y. 

Substitution of the expansions (9) into the equations of motion determines the , 
u ( "\ wf n \ (n > 0), and the y( n ), (n > 0), in terms of the separation profiles x^, u^, and 
w^. However, self-consistency does not pose constraints on the shape of the separation 
profiles themselves. This does not necessarily mean that any separation profile will corre- 
spond to a realistic solution: using the linear heat equation as a model, a velocity profile 
can only correspond to a solution at earlier times if its Fourier transform decays sufficiently 
rapidly with the wave number. For this simple model problem, suitable profiles, which are 
arbitrary close to incorrect solutions, can be found by truncating the Fourier transforms 
of the incorrect profiles at large wave numbers. 

To examine whether modification of the velocity profile in the wall layer is indeed 
possible, the computations were repeated for the case in which the sphere is gradually 
brought to a halt in the time interval 4 < t < 4.5(< ts). The choice t = 4 was made 
because a crude preliminary estimate suggested that ts — 4 was too short a time interval 
for diffusion to reach the particle -qs. The velocity change w r as prescribed as 


W ^ 0 ’^ 1 +exp(277(l - T 2 )) ’ 


(t - 4.25) 
0.25 


(11a, b) 


As can be expected, the results in figure 3 show that the Lagrangian separation profiles are 
dramatically altered near the wall. Similarly figure 4 shows the difference in the Eulerian 
velocity profiles in the wall layer. However, the separation at particle 77 s proceeds exactly as 
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before: the spin of the sphere is brought to a halt but separation continues. At and beyond 
particle S, the Lagrangian velocity profiles with and without spin-down agree within 10 -7 . 

The asymptotic description of the separating layer proceeds in the same manner as 
for the wall layer, with one distinction: the vertical position y is now singular. The reason 
is evident from the integral (2), which turns singular in passing point S. The resolution 
given by Van Dommelen fe Cowley (1989) is to refer y to a reference position y + in the 
separating layer. Many definitions are possible for this reference, but a convenient example 
is: 

w(y + ,t) = 0.05w'(0,t) . (12) 

For this definition y + corresponds to a typical boundary-layer thickness. The continuity 
equation (2) may now be written as: 

, r d v' + 
y = y + /.^i)" y +y ’ 

where y is regular in the separating layer. Thus, when y is replaced by y and v with v = y, 
the description of the separating layer becomes of the same form (9), (10) as for the wall 
layer. 

5. The Emergence Of Logarithmic Terms 

In the previous two sections, the asymptotic expansions for the inner region and the 
wall and separating layers have been found. However, the position ys of the particle 
r]s in the inner expansion, and the reference position y + in the separating layer remain 
undetermined. 

To find ys, the inner expansion (7) can be matched with the wall layer (9), to yield: 

ys(t) ~^ + 2Bln ^ + 0(1) ' (14) 
where the coefficients (3 and B are given in (7ft, c). 

The logarithmic term will lead to a corresponding logarithmic term in the Eulerian 
expansion of the velocity profile (8a), when rewritten in terms of Y: 

G — |(1 — cos/?Y) + - .4 sin /?Y + ... , (15a) 

-2/3B\6t\Un~ + (-2/?Hln(|(l - cos(/?Y)) +4,)|ft|* , (156) 
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where A 0 is a constant depending on the wall layer profile. 

The expansion originally proposed by Hanks Sc Zatnrska (1979) was of the same form 
(15a), but they proposed a different coefficient A, without the logarithmic term: 

Abz ~ 4.6|6t| * . 


To verify that there is in fact a missing logarithmic term, a numerical variable approxi- 
mating A was constructed as: 


_ 7T Y 2 - 3Y, , x 

“ = 2 Y a - Yi ’ (16 “> 

where Yj and Y 2 are the two Y-positions where G = |. For the Banks & Zaturska (1979) 
proposal a — Abz 4- 0(|£t|), while for (156) 


a ~ 




(164) 


According to the present more accurate numerical results, the computed values for a do not 
appear to approach the constant given by Banks Sc Zaturska (1979) (figure 5). If however 
the logarithmic term in (166) is first subtracted from the curve, the required approach to 
a constant does become evident. 

Simpson Stewartson (1982) were the first to point out the presence of logarithmic 
terms based on an Eulerian description of the flow. However, our result B = —0.3957 
agrees poorly with the value —0.457 found by Simpson Sc Stewartson (1982), who obtained 
their value from subtracting two large quantities and fitting of a straight line to the result- 
ing smaller quantity. We submit that our value is independently supported both by the 
apparent convergence of the results in table 3, using the definitions in (76, c), and also by 
its apparent success in eliminating the blow up of the curve figure 5. 

With the inner solution now fully determined by (14), the position of the separating 
layer follows from matching as: 




2tt 

iwT 


+ 0 ( 1 ) • 


(17) 


The 0(1) constant is related to the arbitrariness in the possible definition of the reference 
position y + and the shape of the wall layer profile. 
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6. The First Separated Stages 

Compared to the Eulerian computations, the most remarkable aspect of the La- 
grangian computation is the absence of any numerical difficulties in the integration of 
the momentum equations near the initial separation time ts . The question arises whether 
the Lagrangian solution continues to exist beyond this time. The numerical evidence in 
tables 1 and 3 and figures 6 through 8 indicates that a formal solution does indeed exist 
for a finite range of times tg < t < ty beyond initial separation. Even at 1025 points 
across the boundary layer, there is no sign of convergence difficulties nor of instability. 
And theoretically, at least the coefficient of the second order viscous derivatives in (16) 
and (lc) remains positive, so that the time-like direction remains positive. 

For these reasons were are led to assume that the Lagrangian boundary-layer problem 
has a formal solution beyond initial separation. Whether such a formal solution does have 
a physical meaning is a second and separate question. Certainly at the equatorial plane 
itself, the separation must set up interactive processes which invalidate the further use of 
the original equations (1). 

However, there is a second interpretation for the solution. Van Dommelen (1986) 
computed the Lagrangian solution to the full two-dimensional boundary layer around the 
sphere, a solution which also appears to remain regular at and beyond the separation time 
ts . Now, if it turns out that the interactive effects remain restricted to a small vicinity 
of the equator, say to x = + 0(Re~ a ) with a > 0, then the two half boundary layers 

x < j7r and x > iff would continue to describe the correct asymptotic limit away from 
the equator. These two half boundary layers would near x = iff match with the small 
interactive region. The notion of a limited interaction region seems in line with steady 
descriptions, (Stewartson 1958, Smith & Duck 1977), and the formation of an equatorial 
jet (Dennis &: Ingham 1979, Dennis & Duck 1988). For a limited interaction region, the 
present solution describes the flow in the matching region Re~ Q <C |x — ~n\ >C 1 since it 
is the equatorial limit of the two-dimensional Lagrangian solution. 

Whether or not these arguments apply, the boundary-layer equations are important 
enough by themselves that their possible behaviour is worth study. Note that a solution 
to the boundary-layer equations may be physically relevant in some settings even if it does 
not in apply in other settings (cf. Smith 1982). 

The initial characteristics of the continuity equation beyond time ts were found in part 
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] and are such that the singular behaviour in the continuity equation remains restricted 
to x = Because of the boundary conditions (1 g,j) for x ^ , beyond time t.g the single 
zero for X ( must separate into two zeros, at positions which will be denoted as rjj and 
7 / 2 . Computed values for rjj and 772 are given in table 3 and shown in figure 6. The two 
main boundary layers x < and x > ^7 r each split into three layers near the equatorial 
plane. The layer of particles 77 < 772 remains close to the wall. The middle layer of particles 
Vi < V < V2 penetrates relatively far from the wall on account of the strong growth in y 
near 77] . The upper layer of particles 77 > 772 penetrates still further from the wall because 
of the additional growth in y near V 2 - The middle layer further ejects a finite mass flow 
into the vicinity of the equatorial plane. It may be conjectured that this mass will develop 
into the equatorial jet. 

7. The Settle-Down Of The Lower Stationary Point 

When the time ty beyond the time ts of first separation is approached, new phenomena 
start to show up. The lower stationary point approaches the wall, figure 6, table 3 , and 
the local value of the Lagrangian gradient u^ appears to blow up. Similarly the wall shear 
gradient u^ n rapidly increases, table 1. The “settle-down” of the lower stationary point is 
depicted in figure 6. No significant singular behaviour is evident at the upper stationary 
point, cf. table 3 . 

Since an inviscid flow does not turn singular in the Lagrangian coordinate system, it 
is likely that viscous effects are a primary influence near settle-down. A balance of the 
viscous and convective terms in the Lagrangian boundary-layer equations (1) is consistent 
with an inner coordinate 

E = - — j- , t = t — ty . (18a, 6) 

M 3 

Indeed the U-position Eo of the lower stationary point appears to remain non-zero and 
finite near time ty, figure 7 . Since x^ = 1 at the wall, and vanishes at the stationary 
point, the appropriate asymptotic expansion should be 

l J (E) 

" -j^p • ( 18 c, d) 

This agrees with figure 7 , where the ff-value U 0 at the stationary point and the minimum 
value U m appear to remain finite. In addition the wall shear gradient ratio |u,^|*/|u ) ( 1J |3 
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shown figure 8 appears to remain finite. The correspondingly scaled Eulerian variables are: 

y = |r|*F, u x = -T^r, v=-^ T . (18e, f,g) 

l T l |t|j 

The azimuthal velocity appears to approach a unit value, cf. w >y in table 1 and wj in table 
3. 

Substitution of the inner expansion (18a, b) into the Lagrangian equations (la) yields: 


U + \EU' = x\U" - x,,x' t( U’ + (*'j - x i( x" )(/ , (18ft) 

V = \Ex\ ( , (18*) 

The equivalent Eulerian problem reads 

G" ~(V + \ Y)G' + G 2 -G = 0 , (18;) 

i 

G = V ' , (18fc) 


where accents now denote derivatives with respect to Y. The Eulerian solution is only 
defined below the stationary point E 0 . 

The Runge-Kutta solutions to the inner problems and (18j, k) are shown in 

figure 9. Briefly, the parts 0 < Y < oo and 0 < E < E 0 were found from upward shooting, 
and seeking the least singular solution. The part Eo < E < oo was found from downward 
shooting, starting from a self-consistent asymptotic series truncated at the smallest term. 
The resulting values of the Lagrangian derivatives at Eq are listed in table 4. In the vicinity 
of the stationary point Eo, the solution for x^ may be described to four digits accuracy 
bv a diverging Taylor series of the form 

*,< = F C »( JE - B >) n ■ < 19 ) 

n 

with coefficients given in table 5. 

According to figures 7 and 8, the asymptotic solution of figure 9 is in good agreement 
with the numerical data. In both figures 7 and 8, in order to avoid the difficulty in choosing 
a precise value for time ty , an apparent time difference (t\- — t)* was defined as 

_ 2 

(lv " tr = (few) * ’ G ' (0) = I - 277980908 • ( 2 °) 
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Table 6 shows that this apparent time difference is in good agreement with the real time 
difference near time ty (which also provides more support for the applicability of the Itnnge 
Kutta solution). 

Summarizing the results, it would seem that the Lagrangian solution can be continued 
through separation until the lower stationary point reaches the wall. Beyond that time, 
the limiting singular behaviour of the main two boundary layers x < 5* and x > ^7T near 
the equatorial plane can possibly only be found from a complete integration of these layers. 

8. Concluding Remarks 

Van Dommelen & Cowley (1989) showed that self-consistent unsteady separation pro- 
cesses can be derived bv assuming a smooth Lagrangian solution. To give the strongest 
possible verification that this concept is physically meaningful, in this paper the equatorial 
boundary-layer separation for the impulsively spun sphere was recomputed. With only one 
spatial dimension, this case allows excellent numerical accuracy. Our scheme appears the 
most accurate yet; its separation time is accurate to seven digits before using Richardson 
extrapolation. 

Even when we used over a thousand points across the boundary layer, we could not 
observe any deviations from the smooth Lagrangian solution proposed by Van Dommelen 
cV Cowley (1989). Derivatives up to third order could easily be determined to five digits 
accuracy, cf. table 3. Derivatives of still higher order would be difficult to evaluate, 
but they play no significant part in the final separation structure. Moreover, singular 
behaviour of the higher order derivatives would tend to render the evaluation of the lower 
order derivatives more difficult, and we observed no evidence of that. 

Numerical continuation of the boundary-layer solution beyond the time of first sepa- 
ration showed that the wall vorticity layer disappears in a finite time. Whether a similar 
process occurs for the non-interactive solution of asymmetric two- or three-dimensional 
separation, in which the separation structure is in motion compared to the wall, remains 
unknown. 

During parts of this investigation, the author was supported by the ONR, the AFOSR, 
and ICOMP, NASA Lewis. 
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Table 1. The wall shear at the equatorial plane of the spinning sphere. The top tabulation 
yields values of the wall shear w >y and the bottom tabulation values for the wall 
shear gradient u xy . 


mesh 

points 

time 

step 

ts : 

tv : 

: 129 

0.025 

4.575676 

5.5418 

257 

0.0125 

4.575643 

5.5461 

513 

0.00625 

4.575634 

5.5473 

1025 

0.003125 

4.575632 

5.5476 

Banks 

Sz 

Zaturska 

4.5758 

Simpson 

Sz 

Stewartson 

4.57446 

t = 0 

-0.56415 

-0.56418 

-0.56419 

-0.56419 

-0.564190 


0.25 




-1.12666 

-1.12668 


1 


-0.55025 


-0.55025 

-0.55029 

-0.5502 

2 


-0.35771 

-0.35771 

-0.35771 

-0.35773 

-0.3577 

3 


-0.24414 

-0.24416 

-0.24417 

-0.24418 

-0.2441 

4 

-0.14417 

-0.14445 

-0.14452 

-0.14453 

-0.1446 

-0.1445 

4.3 

-0.11450 

-0.11487 

-0.11496 

-0.11498 

-0.1150 

-0.1149 

4.4 

-0.10469 

-0.10508 

-0.10518 

-0.10521 


-0.1052 

4.5 

-0.09494 

-0.09537 

-0.09547 

-0.09550 

-0.0955 

-0.0954 

4.52 

-0.09300 

-0.09343 

-0.09354 

-0.09357 

-0.0936 

-0.0935 

4.54 

-0.09106 

-0.09150 

-0.09161 

-0.09164 

-0.0917 

-0.0916 

4.56 

-0.08913 

-0.08958 

-0.08969 

-0.08971 

-0.0897 


ts 

-0.08761 

-0.08807 

-0.08819 

-0.08821 



5. 

-0.04839 

-0.04859 

-0.04866 

-0.04867 



5.25 

-0.02632 

-0.02653 

-0.02661 

-0.02663 



5.39 

-0.01424 

-0.01441 

-0.01449 

-0.01452 



5.47 

-0.00741 

-0.00736 

-0.00745 

-0.00748 



5.5 

-0.00527 

-0.00463 

-0.00472 

-0.00475 



5.537 



-0.00115 

-0.00118 



t = 0 

-0.40990 

-0.41000 

-0.41002 

-0.41003 

-0.41003 


0.25 




-0.20527 

-0.20527 


1 


-0.41855 


-0.41858 

-0.41855 

-0.4186 

2 


-0.63387 

-0.63396 

-0.63398 

-0.63395 

-0.6340 

3 


-0.89600 

-0.89618 

-0.89622 

-0.89619 

-0.8963 

4 

-1.40385 

-1.40486 

-1.40513 

-1.40520 


-1.4058 

4.3 

-1.71445 

-1.71522 

-1.71543 

-1.71549 


-1.7166 

4.4 

-1.85470 

-1.85794 

-1.85810 

-1.85814 


-1.8594 

4.5 

-2.03105 

-2.03121 

-2.03127 

-2.03128 


-2.0330 

4.52 

-2.07037 

-2.07043 

-2.07047 

-2.07048 


-2.0722 

4.54 

-2.11148 

-2.11143 

-2.11143 

-2.11144 


-2.1133 

4.56 

-2.15450 

-2.15432 

-2.15429 

-2.15429 



ts 

-2.18962 

-2.18927 

-2.18920 

-2.18918 



5 

-4.132 

-4.110 

-4.104 

-4.103 



5.25 

-9.026 

-8.903 

-8.873 

-8.865 



5.39 

-22.34 

-21.63 

-21.45 

-21.41 



5.5 

-151 

-129 

-125 

-124 



5.537 



-1231 

-1171 
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Table 2. The quantity |u >x | m a t — (ts — t) 1 according to the present four meshes and 
according to Simpson Sc Stewartson (1982). 



129 

X 

0.025 

257 

X 

0.0125 

513 

X 

0.00625 

1025 

X 

0.003125 

Simpson 

& 

Stewartson 

t = 4 

0.035150 

0.035127 

0.035121 

0.035120 

0.0352 

4.3 

0.017546 

0.017511 

0.017501 

0.017499 

0.0176 

4.4 

0.011318 


0.011243 

0.011240 

0.0112 

4.5 

0.00504 

0.00491 

0.00488 

0.00486 

0.0045 

4.52 

0.00381 

: 

0.00362 

0.00359 

0.00359 

0.0027 



Table 3. Stationary point trajectory figure 6 and corresponding Lagrangian quantities for 
the four increasingly finer meshes of table 1. Note that these values do not 
Richardson extrapolate; they were found from second order interpolation between 
the mesh points. However, at the time ts of first separation, the stationary point 
happens to fall nearly exactly on a mesh point of the finest two meshes. The four 
values for each quantity are in order of increasing accuracy. 


t 

V 

w 

u ,« 

x ,tw 

u ,«n 

x ,(vvv 

ts 

0.97189 

0.672593 

-0.369784 

0.689301 

-0.04000 

-1.1328 


0.97185 

0.669262 

-0.369789 

0.689116 

-0.04019 

-1.1291 


0.97186 

0.669249 

-0.369789 

0.689064 

-0.04026 

-1.1276 


0.97188 

0.669240 

-0.369790 

0.689062 

-0.04026 

-1.1274 

t 

Vi 

Wi 

u .<i 

V2 

W 2 

CN 

3 

5 

0.41478 

0.92405 

-0.46260 

1.81402 

0.39848 

-0.40755 


0.41507 

0.92389 

-0.46196 

1.81422 

0.39838 

-0.40763 


0.41513 

0.92385 

-0.46180 

1.81426 

0.39837 

-0.40764 


0.41515 

0.92385 

-0.46176 

1.81427 

0.39837 

-0.40764 

5.25 

0.28209 

0.97193 

-0.7501 

2.09452 

0.33781 

-0.41589 


0.28281 

0.97167 

-0.7454 

2.09465 

0.33776 

-0.41597 


0.28300 

0.97159 

-0.7442 

2.09468 

0.33774 

-0.41599 


0.28304 

0.97158 

-0.7439 

2.09469 

0.33774 

-0.41600 

5.39 

0.20346 

0.98922 

-1.394 

2.23830 

0.31096 

-0.41926 


0.20473 

0.98900 

-1.370 

2.23844 

0.31089 

-0.41937 


0.20494 

0.98900 

-1.364 

2.23846 

0.31088 

-0.41939 


0.20500 

0.98900 

-1.363 

2.23846 

0.31088 

-0.41939 

5.47 

0.1435 

0.99643 

-3.025 

2.31751 

0.29725 

-0.42090 


0.1463 

0.99610 

-2.874 

2.31763 

0.29720 

-0.42101 


0.1469 

0.99606 

-2.841 

2.31765 

0.29719 

-0.42103 


0.1470 

0.99605 

-2.834 

2.31766 

0.29718 

-0.42103 

5.5 

0.1124 

0.99831 

-5.33 

2.34681 

0.29231 

-0.42152 


0.1161 

0.99811 

-4.89 

2.34686 

0.29231 

-0.42158 


0.1171 

0.99805 

-4.78 

2.34686 

0.29231 

-0.42160 


0.1173 

0.99804 

-4.75 

2.34687 

0.29231 

-0.42161 

5.537 

— 

— 

— 

— 

— 

— 


0.0573 

0.99978 

-24.0 

2.38259 

0.28648 

-0.42227 


0.0581 

0.99977 

-23.2 

2.38259 

0.28648 

-0.42228 
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Table 4. First few Lagrangian derivatives of the solution figure 9 at the stationary point 
Ett. 


solution 

x ,e 

x !e 

x" 

x ,« 

U 

U' 

V" 

E < E o 

0 

0.92099 

2.015 

-0.27450 

0.13975 

1.02 

E > Eq 

0 

0.92091 

2.02 

-0.27450 

0.13977 

1.02 


Table 5. First 30 coefficients C n in the self-consistent Taylor series expansion (19) to the 
solution figure 9. They show no evidence of a non-zero radius of convergence. 
The tabulated results do not depend critically on the precise values of E 0 and 

c,. 


O.OOeOO 

— 0.92e00 

O.lOeOl 

— 0.56e00 

— 0.23e00 

0.07e00 

0.65e00 

0.24e01 

0.11e02 

0.73e02 

0.65e03 

0.74e04 

0.10e06 

0.17e07 

0.33e08 

0.72e09 

O.lSell 

0.48el2 

0.14el4 

0.47el5 

0.16el7 

0.63el8 

0.26e20 

0.11e22 

0.52e23 

0.25e25 

0.13e27 

0.72e28 

0.41e30 

0.25e32 

0.16e34 



Table 6. Comparison between the apparent time from settle-down as defined in (20) and 
the true value, (mesh 1025 x 0.003125; ty = 5.54760) 


v/ty — t 

0.7218 

0.5402 

0.3953 

0.2781 

0.2180 

0.1661 

0.1503 

0.1326 

0.1166 

0.1030 

0.0748 


y/ty ~ t 

0.7400 

0.5455 

0.3970 

0.2786 

0.2182 

0.1661 

0.1503 

0.1327 

0.1166 

0.1030 

0.0748 
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FIGURE 1. - LAGRANGIAN BOUNDARY-LAYER PROFILES IN THE EQUA- 
TORIAL PLANE OF THE IMPULSIVELY SPUN SPHERE AT THE TIME 
t s = 4,575632 THAT SEPARATION STARTS. HERE w IS THE AZI- 
MUTHAL VELOCITY COMPONENT AND u £ THE LAGRANGIAN GRADIENT 
OF THE MERIDIONAL VELOCITY U. FURTHER x, IS THE LAGRANGIAN 
GRADIENT OF THE POLAR PARTICLE POSITION x, THE SEPARATION 
PARTICLE n s = 0.97188 IS INDICATED BY A SOLID DOT. THE 
CONDITION OF VANISHING x > IMPLIES A SINGULAR EULER I AN VE- 
LOCITY PROFILE. 





FIGURE 2. - EULER I AN VELOCITY PROFILES WHEN THE TIME OF FIRST 
SEPARATION IS APPROACHED. Y IS A SCALED DISTANCE (3b) FROM 
THE WALL AND G IS THE SCALED EULER IAN VELOCITY GRADIENT u x 
DEFINED IN (43). NOTE THAT NEAR SEPARATION MOST OF THE 
EULER IAN VELOCITY PROFILE CORRESPONDS TO ONLY A SMALL VI- 
CINITY OF THE STATIONARY POINT r) s OF THE LAGRANGIAN PROFILE 
FIGURE 1. 
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FIGURE 3. - LAGRANGIAN PROFILES WHEN THE SPIN OF THE SPHERE 
IS SMOOTHLY BROUGHT TO A HALT. COMPARISON WITH THE PRE- 
VIOUS LAGRANGIAN SEPARATION PROFILES FIGURE 1 SHOWS THAT 
THE VELOCITY CHANGE GENERATES A SUB-LAYER AT THE WALL. 
HOWEVER, THIS LAYER DOES NOT REACH THE PARTICLE n s IN TIME 
TO HALT THE SEPARATION PROCESS. 




FIGURE - EULER IAN PROFILES CORRESPONDING TO FIGURE 3 



FIGURE 5. - PRESENCE OF HIGHER ORDER LOGARITHMIC TERMS 
IN THE EXPANSIONS FOR THE IMPULSIVELY SPUN SPHERE. 
ACCORDING TO THE ORIGINAL PROPOSAL OF BANKS AND ZATURSKA 
(1979), THE UPPER CURVE SHOULD LINEARLY APPROACH THE 
VALUE 4.6 INDICATED BY THE DOT. BUT ACCORDING TO THE PRE- 
SENT RESULTS THE UPPER CURVE CONTAINS A LOGARITHMIC TERM, 
AND ONLY AFTER SUBTRACTION OF THIS TERM THE LOWER CURVE 
APPROACHES A CONSTANT VALUE. SEE TEXT FOR A COMPARISON 
WITH THE WORK OF SIMPSON AND STEWARTSON (1982). 



FIGURE 6. - TRAJECTORY OF VANISHING X,£ FOR THE IMPULSIVELY 
SPUN SPHERE IN THE LAGRANGIAN EQUATORIAL r). t-PLANE. IN 
CONTRAST TO THE EULER I AN CASE, THE LAGRANGIAN SOLUTION DOES 
NOT APPEAR TO TERMINATE AT THE TIME tg THAT SEPARATION STARTS. 
BUT NEAR THE SETTLE-DOWN TIME t v = 5.5476 WHEN THE LOWER STA- 
TIONARY POINT ATTACHES ITSELF TO THE WALL, THE LAGRANGIAN SOL- 
UTION TURNS SINGULAR ALSO. 
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FIGURE 7. - BEHAVIOR OF THE LAGRANGIAN SOLUTION WHEN THE SETTLE- FIGURE 8. - THE AGREEMENT IN WALL SHEAR WITH THE ASYMPTOTIC 

DOWN TIME t v IS APPROACHED. THE SCALED VARIABLES ARE DEFINED SOLUTION FIGURE 9 WHEN THE SETTLE-TIME t y IS APPROACHED. 

IN THE TEXT. THE SOLID DOTS DENOTE THE VALUES 0.5961, 0.2824, THE SOLID DOT DENOTES THE VALUE 1.101615 ACCORDING TO THE 

AND 0.2745 ACCORDING TO THE ASYMPTOTIC SOLUTION FIGURE 9. ASYMPTOTIC SOLUTION FIGURE 9. 



0 1 0 1 


FIGURE 9. - PROPOSED ASYMPTOTIC SOLUTION IN THE VICINITY OF SETTLE-DOWN OF THE LOWER 
STATIONARY POINT IN FIGURE 6. THE EULER I AN SOLUTION 0 — Y < 00 IS EQUIVALENT TO THE 
LOWER LAGRANGIAN SOLUTION 0< E < E n . 
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